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Abstract 

We prove some sharp Hardy type inequalities related to the Dirac operator by 
elementary, direct methods. Some of these inequalities have been obtained previ- 
ously using spectral information about the Dirac-Coulomb operator. Our results are 
stated under optimal conditions on the asymptotics of the potentials near zero and 
near infinity. 
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1 Introduction and main results 



The uncertainty principle is without any doubt a fundamental attribute of 
quantum mechanics [17]. In the case of the Laplacian it states that for all 

* © 2003 by the authors. This paper may be reproduced, in its entirety, for non- 
commercial purposes. 

1 PICS (CNRS, Paris-Bilbao) and European Programs HPRN-CT # 2002-00277 
& 00282. 

2 Work partially supported by U.S. National Science Foundation grant DMS 00- 
70589. 



7 February 2008 



functions / G C^°(M 3 ), 

This inequality is also known as Hardy's inequality. By scaling, the power of 
the potential is seen to be optimal but also the constant | cannot be improved. 
However, it is still possible to improve the inequality by adding lower order 
terms. 

In recent years there has been a great effort to find optimal improved Hardy- 
type inequalities in the case of the Laplacian. The pioneering work in this 
direction is due to Brezis and Vazquez [7] in the case of Dirichlet boundary 
conditions, and to Lieb and Yau [16] in the case without boundary conditions. 
Further improvements have been obtained in [1,2,3,4,5,6,11,18]. 

An analogue of this inequality for a relativistic version of the Schrodinger equa- 
tion where the Laplacian is replaced by \J — A is an inequality due to Kato : 

(f,V^Af)>- [J-^-dx (2) 

v ' 7T JR 3 \X\ 

(see [14,15]). In this inequality the power and constant are again optimal. 
An immediate consequence of this inequality is that the relativistic model of 
the hydrogenic atom with kinetic energy \J — A and with nuclear charge Z is 
stable if and only if v := Za < ir/2, where a ~ 1/137. . . is the fine structure 
constant. Lieb and Yau [16] also discovered some generalizations for \J — A to 
balls. 

The Dirac relativistic hydrogenic atom is stable only if v < 1. The Dirac 
Hamiltonian is unbounded from below and instability has to be interpreted 
in the different, more subtle, sense of a breakdown of selfadjointness of the 
Dirac-Coulomb Hamiltonian. If the Coulomb singularity is smeared out, then 
the threshold for stability is reached when the lowest eigenvalue in the gap 
reaches the upper bound of the negative continuum. This happens in general 
for larger values of v. In the case of the Dirac-Coulomb Hamiltonian, the 
stability is a consequence of the following Hardy-type inequality. 

Theorem 1 [9] Let <r = (o"j)j = i i2 ,3 be the Pauli-matrices : 

*1 = (2 I), *2 = (S o*), <73= (o 

Then for every up G H\R 3 , C 2 ), 
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As in (1) and (2), the powers of \x\ and the constants are optimal. Inequal- 
ity (3) has been established using a characterization of the eigenvalues of a 
self-adjoint operator in a gap of its essential spectrum by means of a particular 
min-max. See below and refer to [9] for more details. For other results on min- 
max characterizations of the eigenvalues of Dirac operators, see [10,13,12,8]. 

By scaling, if we replace </?(•) by £~ 1 ip{e^ 1 -) and take the limit e — > 0, (3) 
implies that for all tp G H 1 ^ 3 , C 2 ), 

/ Id \a ■ Vip\ 2 dx > [ ^-dx. (4) 
it 3 ' 1 ii 3 \x\ 

This inequality is slightly generalized form (take (p = (g, 0) and consider in- 
dependently the cases where g takes either real or purely imaginary values) of 
the following inequality : for all g G if 1 (M 3 , C), 

/ \x\ \'Vg\ 2 dx> [ dx , (5) 

jR 3 JR3 \x\ 



which is itself equivalent to (1) : take / = y\x\g. Note that the largest 
space in which Inequality (3) holds is larger than i? 1 (R 3 ,C 2 ) and contained 
in H^ 2 (R 3 ,C 2 ). For more details, see [9]. 



In [9], the proof of (3) has been carried out by using explicit knowledge on 
the point-spectrum of the Dirac-Coulomb operator H u :— — % oc • V + /3 — ^ , 
where the matrices (3, G A^4 X 4(C), k = 1, 2, 3, are defined by 

«*=U ^=(o -id)- 

Id is the identity matrix in C 2 and v is a real parameter taking its values in 
(0, 1). It is well-known [19] that for any v G (0, 1) H u can be defined as a self- 
adjoint operator with domain V v satisfying : /^(R 3 , C 4 ) C V u C # 1/2 (R 3 , C 4 ) 
and spectrum 

a(H u ) = <j ess (H ) U |Ai, A2, • • • I , a ess (i/ ) = (-oo,-l]U[l,+oo) , 

where {A£}fc>i is the nondecreasing sequence of eigenvalues of H u , all contained 
in the interval (0, 1) and such that : 

Ai = Vl - ^ 2 , lim A^ = 1 for every 1/ G (0, 1) . 

For a large set of potentials V with singularities not stronger than Coulombic 
ones, more precisely, for all those satisfying : 

lim V(x) =0 and - ^ - ci<V <c 2 = sup(V) , (6) 

|a;|^+oo \x\ 
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with v G (0, 1), Ci, c 2 G R, Ci, c 2 > 0, Ci + c 2 — 1 < \A — z/2 > the following 
result was proved in [9] : 

Theorem 2 [9] Let V be a radially symmetric function satisfying (6) and 
define Ai(V) as the smallest eigenvalue of H + V in the interval (—1,1). 
Then for all ip G H^R^C 2 ), 

Inequality (7) is achieved by the large component, i.e., the two-spinor made 
of the first two complex valued components of the four-spinor, of any eigen- 
function associated with Ai(V). In particular if V — —j^, v G (0, 1), we get 

Corollary 3 [9] For any v G (0, 1), for all up G H\R 3 , C 2 ), 

r ( + (i-vt^) \<pA d X > v [ ^ dx . (8) 

This inequality is achieved in L 2 (R 3 , \x\~ l dx) 4 . Inequality (3) is obtained 
from (8) by taking the limit v — > 1. Theorem 1 is therefore a straightforward 
consequence of Corollary 3. Note that (3) is not achieved in L? (R 3 , |x| -1 ofo;) . 

The aim of this paper is twofold. On the one hand, we give a direct analytical 
proof of Theorem 1 which does not use any a priori spectral knowledge on 
the operator H u . On the other hand, we prove more general inequalities by 
showing that for some continuous functions W > 1 a.e. and constants R > 
and C(R) < 0, the inequality 

holds for all (p G i? 1 (R 3 ,C 2 ). Here /zr is the surface measure induced by 
Lebesgue's measure on the sphere Sr := {x G R 3 : \x\ = R}. Note that this 
inequality is relevant for the Dirac operator with potential V — rjr. Improved 
inequalities like (3) or (9) with the operator <r • V replaced by V can easily 
be obtained by considering separately the real and the imaginary parts of the 
components of the two-spinors. 

We are actually interested in understanding for which functions W Inequal- 
ity (9) holds and what is the optimal behavior of the function W near or 
near +oo. By optimal at s = or s = +oo, we mean optimal at each order in 
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the sense that we look for a expansion of the form W = 1 + J2t=i °k Wfc such 
that at each order ko, < W ko +i = o(W — 1 — J2 k =\ °k W k ) and c ko+ i is the 
largest possible constant. What we are going to prove is a little bit involved: 
it is not clear that c k is a constant (see Appendix C) and we are only going to 
prove that the maximum of its lim inf is achieved. As in Hardy-like inequalities 
for the Laplacian operator on balls centered at the origin, we will see that the 
optimal behavior of the function W near is a logarithmic perturbation of 
the constant 1. More precisely the optimal behavior near the origin is given 
by functions of the form 

1 oo 

WM = 1 + -£^(M) 2 ---A^(kl) 2 , 

y k=l 

where Xi(s) := (a — log(s)) -1 for some a > 1, X k (s) := Xi o X k -i). The 
functions X k and W°° are well defined for \x\ — s < e a ~ l (see Appendix A for 
basic properties of the functions X k ). These asymptotics are optimal in the 
above sense, with W k = Xx(|a;|) 2 • • -X fc (|a;|) 2 and liminf c k = |. On the other 
hand, as \x\ — > +oo, the optimal growth for W is given by \x\, i.e., the first 
term in the l.h.s. of (9) does not help. 

Theorem 4 Assume that for some R > Inequality (9) holds for every spinor 
(p G Cq°(R 3 ,C 2 ) ; where W is a radially symmetric continuous function from 
]R + to IR + . Assume moreover that W(0) > and W is nondecreasing in a 
neighbourhood of + . Then W(0) < I, 

lim sup W(s)/s < 1 

and for all k > 1, 

liminf (W(s) -1-1 j^Xfis) ■ • • X 2 ( S ))xr 2 ( S ) • ■ ■ X^(s) < i . (10) 

As soon as W ^ 1, C{R) must be negative. Moreover, there are continuous 
functions W > 1 for which (9) holds with some R > 0, C(R) < 0, such that 
lim s ^ +00 W(s)/s = 1 and (10) holds with equality for all k > 1. 

Note that this result is independent of the particular value of a > 1 which ap- 
pears in the definition of the functions Xj,. The fact that we have to introduce 
a discontinuity at some s = R is not contradictory with the known results for 
the usual Hardy inequality, for which bounded domains are considered and q 
is taken large enough. 

In Section 2 we give a direct analytical proof of Corollary 3 together with 
several auxiliary results. We recall that Theorem 1 is a straightforward con- 
sequence of Corollary 3. Section 3 is devoted to the proof of Theorem 4. 
Properties of the functions X k , an existence result for a singular ODE needed 
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for Theorem 4 and an example illustrating why we have to consider a liminf 
in this theorem are given in three appendices. 



2 Proof of Theorem 1 

In this Section, we actually prove Corollary 3, which is slightly more than 
Theorem 1. 

First we fix some notations. The spinor (p = ( %\ ) takes its values in C 2 and 
by |(/?| 2 , |V(^| 2 and |cr ■ V(/?| 2 we denote, respectively, the quantities |y?i| 2 + 
I V^2 1 2 , S| =1 (|9 fc v9i| 2 + |<9 fc y? 2 | 2 ) and |<9 3 y?i + <9i</? 2 — i <9 2 <^ 2 | 2 + |#i<Pi + * - 
d 3 ip 2 \ 2 - Further, we notice that the Pauli matrices are Hermitian and satisfy 
the following properties : 

Uj a k + a k Oj = 2 5 jk Id , V j, k — 1, 2, 3 . 

With a standard abuse of notations, each time a scalar 5 appears in an identity 
involving operators acting on two-spinors, it has to be understood as S Id, 
where Id is the identity operator. 

On the other hand, for all vectors a, b G C 3 , we have 

(cr • a) (cr • b) = a ■ b + i cr ■ (a x b) . 

Applying this formula to a = x and b = iV, we obtain the following expres- 
sion for the commutator of cr ■ V and ax: 

[cr • V, cr ■ x] = V-x-x-V + 2cr-L = 3 + 2cr-L, 

where L = — i x A V is the orbital angular momentum operator. The main 
point to note here is that L acts only on the angular variables. 

For simplicity, for any function h : R + — > R, we denote the functions x \— > 
h{\x\) and x hh> by h and h! respectively. Now if such a function h is 

differentiable a.e. in IR + and continuous in [0, R) U (R, +oo), we have 

[o- • V, (cr-x)h}= \x\ ti + (3 + 2 a ■ L) h + R [h] R 5 R 

= 2 (1 + cr ■ L) h + h + \x\ ti + R [h] R 5 R , 

where by [ti\ R := h(R + ) — h{R ) we denote the possible jump of h at R and 
5 R is the Dirac delta function at r = R, in spherical coordinates. 

The spectrum of the operator 1 + cr • L is the discrete set {±1,±2, •••} 
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(see [19]). This can be seen by noticing that 

l + a-L = J 2 -L 2 + -, J = L + . 

4 2 

Then, the fact that the spectrum of J 2 (resp. L 2 ) is the set + 1) ; j = 
|, §, . . . } (resp. {£(£ + 1) ; £ = j ± |, j = |, |, . . . } proves the above result. 
The main point here is that is not in the spectrum of 1 + a ■ L. If we denote 
by X + (resp. X_) the positive (resp. negative) spectral space of 1 + er ■ L, and 
by P± = \{l ± the corresponding projectors on H 1 (M 3 ,C 2 ), for all 

up e H^R 3 , C 2 ), for all h as above, 

(<p+, [tr ■ V, (a- ■ x) h\ <p+) > [ (3h + \x\ h')\ v+ \ 2 dx + R [h\ R [ \<p + \ 2 dpi R , 

JR 3 JS R 

(</?_, [<r • V, (a ■ x) h] </?_) < / |a;| ti)\(p_\ 2 dx + R[h] R \ip_\ 2 dfi R , 

JR 3 ' Vs fl 

where ip± := P±<£. 

By Cauchy-Schwartz' inequality, for any measurable function g : IR + — > M + , 
/ (3/i+ |x| /?/W+| 2 Gfe< / 9 \(tr ■ V)if + \ 2 dx 

(ii) 

\ x\ 2 h 2 

9 ' ' Js R 



f \x\ 2 h 2 f 

+ / \ip + \ 2 dx - R[h} R \(p + \ 2 d[i R , 

Ji 3 a Js R 



where again we abbreviate g(\x\) by g. 
Define now W and m : M + — > M by 



and assume that W is positive on IR . With the same notation as above, we 
can rewrite (11) as 

f 1 

j R3 -(2m + sm'-(s + W)m 2 )^ xi \ V+ \ 2 dx 

(12) 

\x\ 



Similarly, for <y?_ := P-ip, we find 
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c 1 

/ - f 2m - s m - (s + W) m 2 ) \ip_\ 2 dx 

(13) 

< / l( °|' V) ^ 1 dx + R[h] R [ \<p-\ 2 dpt R . 

Note that for any measurable radial function 6, the spaces X + and X_ are 
also orthogonal in L 2 (R 3 ; b(\x\) dx). Moreover, we have 

Lemma 5 P_(cr- vfp + = P+(o--V) 2 P_ = in PT^R 3 , C 2 ). 

Proof. A direct computation shows that the anti-commutator {cr • V, 1 + cr ■ 
L} = 0, i.e., cr • V anticommutes with 1 + cr • L. Hence (cr • V) 2 commutes 
with 1 + a ■ L. 

Now, let $ ± G X ± be two eigenf unctions of 1 + cr • L with eigenvalues A ± , 
A_ < < A+. Then, 

((cr ■ V)$_, (cr ■ V)$+) = ^ ($_, (cr • V) 2 (l + cr ■ L)$ + ) 

= ^($-,(l+cr- J L)(cr.V) 2 $ + ) 
= ^((l+cr. J L)$_,(cr.V) 2 $ + ) 
= ^-($_,(cr.V) 2 $ + ) , 

which is impossible except if ((cr • V)$_, (cr • V)$+) = 0. □ 

Adding (12) and (13), we get the following result. 

Proposition 6 Let W be a positive measurable function on R + and con- 
sider two functions m± : R + — > R sttc/i £nai £ne maps s \— > m±{s)/s are 
continuous on [0,R) U (P, +oo) and differentiate a.e. on R + . Tnen /or any 
V? G P^C 2 ), 



(g-v)y| 

1 + ^ 

1 + 1*1 



± 



Sr 



> E / R3 \ ( 2 ™± ± " (* + m|) |s=|x| |P ± ^| 2 ^ . 
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In order to prove Inequality (9), we have to find two functions, m + and m_, 
and a continuous function W > 1 such that for s G M + a.e., 

2m± ± s m'± - (W + s) m\ > W - s . 

This means 



i < w < ta^±i^± 2 «"'±t» (14) 



2m± ± s m'_j_ — s m| + s 
1 + m± 

Moreover W has to be as large as possible (near the origin and near infin- 
ity), i.e. optimal in the sense of Section 1. Then (9) follows with C(R) = 
-max([m_] R , [m + \ R ). 

In the sequel, for every function m as in Proposition 6, we will use the notation 

W ±,m := 2m±Sm/ ~! m2 + S . (15) 
1 + m 2 K ' 



Proof of Theorem 1 and Corollary 3. This is simply done by choosing m + = 
m_ = 1 or m + = m_ = 1 -^^- v ' i i n Proposition 6. Since in both cases the 
functions s — > m± / s are continous in M + , there is no surface integral term in 
those inequalities : C(R) = 0. □ 



Remark The above arguments leading to the proof of Theorem 1 can also 
be viewed as a completing the square strategy : for all if G // 1 (1R 3 ) ; h, g G 
C(M + ,lR + ) ; h differentiate a.e. in M, + , it is clear that 



I 

Jr- 



(a ■ V)P ±V ? ± 



(cr • x)h 

V9 



dx > . 



Expanding the squares in the above expressions, integrating by parts the cross 
terms and adding the two inequalities, we find (9) ifW satisfies (14)- 



3 Proof of Theorem 4 



3.1 Direct estimates 



We start this section by finding the optimal behavior near the origin and near 
infinity for continuous functions W for which (9) holds, where W is given 
or not by some function m as defined in (15). This is done in a series of 
intermediate results. 
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Lemma 7 Let W be any function satisfying W > 1 on M + and for which 
Inequality (9) holds. Then, necessarily, 

Wis) 

lim W(s) = 1 and limsup < 1 . 

s->0+ s-»+oo S 



Proof. By assumption, W > 1. If we had VT(0) > 1, it would be easy to 
contradict the fact that 1 is the best constant in (4) and (5). As s — > +oo, 
the result follows from the simple observation that by scaling we can easily 
construct functions (p n := n~ 3 / 2 (p(- /n) such that 

/ \(f n \ 2 dx = 1, / \cr • Vy? n | 2 dx — > as n — > +oo , 

JR 3 JR 3 

so that the gradient term does not play any role. □ 

Proposition 8 Let m± £ C([0,5)) for some 5 > 0. Consider W ±,m± defined 
according to (15) and let W := min(W + ' m+ , W~' m -) be a function for which 
W > 1. Then (9) holds, m±(0) = 1, m± > 1 in a neighbourhood of s = + 
and 

limsup (m±(s) — 1) log(s) < +oo . (16) 

Proof. We prove this for m + , the proof for m_ being identical. Let us write 
m + — n + 1. Then (14) is equivalent to 

< ^_ 1= ^-" 2 -2 a n- 3 n 2 

2 + 2n + n 2 v y 

From this inequality, we infer that n' > n(n + 2) + n 2 /s, so that there 
are two possibilities for the behavior of n near : either n is monotone 
and lim s _ >0 + n(s) = a G (— oo,+oo], or n oscillates near in the inter- 
val (—2,0). The latter case is impossible because on the sequence of local 
minima approaching 0, the r.h.s. of (17) would eventually be negative. So 
lim s ^ + n { s ) — a £ {— °°) +oo] and if a ^ 0, for s > small, 

s n ~ n 2 , 

which by integration implies that near 0, 

for some constant C > 0, a contradiction. Hence, necessarily, a = and the 
result follows from (18), which still holds true when a = 0. □ 

Next we prove the following asymptotic result : 
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Lemma 9 Let A denote the class of the functions n, continuous in the inter- 
val [0, 5) for some 5 > 0, and such that n(0) = 0. Then, for all k > 1, 

sup jliminf (W( S )-n 2 ( S )-i ^Xfo) ■ • • X 2 ( S ))x r 2 ( S ) • • -X^s) J = 1 . 



The fact that we are dealing with a liminf and not a limsup may look surpris- 
ing at first sight. However, an upper limit cannot be expected, as it is shown 
by the example given at the end of the paper. 

Proof. For < s < 1 and Xi(s) := (a — logs) -1 , define implicitely m by 

1 / ( 1 \\ 

n(s) = -Xi(s) 1 -2m 



2 V \Xi(s)JJ ' 

Then 

.. . 9 . . 1 tn'At) — n 2 (t) 1 . 

, rfM -„. ( .) = _ + ^ t 6(li+oo) . 



Next, for all k > 1 and s > 1, let us define again n fc+1 in terms of n k by 
:= Y t { 2n ^ " X ) ' * = xilJT) G (1 ' +00) • 



Then 



At 2 t 2 

Hence, for every k > 1 and every < s < 1, with z = 1/X k (s), we have 
S n'( S )-n 2 ( S ) = \j2Ms) 2 ■ ■■X j (s) 2 + Ms) 2 • • ■ X k (s) 2 (zn' k (z) - n 2 k {z)) . 



1) Choosing n k = delivers a function n(s) with 

4 3=1 

Note that in this case, n(s) = X^Li Ai(s) • • ' Xj(s) ( see Appendix A for more 
details). This shows that 



sup 

A 



jliminf (sn'(s)-n 2 (s)-\ ^X^s) • • • X 2 ( S ))x r 2 ( S ) ■ ■ -X fe - 2 ( S )| > 1 . 

(19) 
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2) Let now n be any function in A. For every k, 

liminf (tn' k — n k ) < . 

If the above limit was to be larger than 0, say some constant b > 0, then, 
integrating the inequality tn' k > n\ would show that n k tends to at infin- 
ity while on the other hand, integrating tn' k > 6/2 would show that n k is 
unbounded near infinity, which provides an obvious contradiction. □ 

Corollary 10 Let W be as in Proposition 8. Then (10) holds and the optimal 
asymptotic behavior near the origin is achieved. 

Proof. Close to s = + , the fact that 

s n' — n 2 — 2 s n — sn 2 1, . 2x 

~ - (sn — n ) 

2 + 2n + n 2 2 v ' 

immediately provides 



liminf (W{s) j^Xlis) ■ ■■X 2 {s))x^{s) • • ■ X l 



k+l\ S ) < Q 



and the optimal behaviour is achieved, for instance, by 

1 +oo 

W = W :=mm(W + ' 1+n ,W-' 1 - n ) , n := T^s)-^) , 

2 3=1 

so that 

1 +oo 

W(s) = l + -Y / X 2 (s)---X 2 (s) + o(s) 
(see Appendix A for more details). □ 



3.2 Estimates based on improved Hardy inequalities for the Laplacian 



The above arguments show that the optimal growth near and near infinity 
for any function W generated (as above) by functions m±, continuous near the 
origin and near infinity, and for which (9) holds, is given by (10) with equality 
for each k > 1, as in the statement of Theorem 4. On the other hand, the 
optimality near infinity was established in Lemma 7. However, it remains to 
prove that there is no function W - not given by (15) - with higher growth at 
the origin. This amounts to prove that there is no radial function W with more 
singular asymptotics near the origin and for which the differential problem 

sm = sm 2 — s — 2m + (1 + m 2 ) W , m(0) = 1 , 
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cannot be solved for some function m, continuous at 0. The rest of this section 
is devoted to this question. 

Step 1 : We first remark that in this problem the angular variables do not play 
any role : only radially symmetric spinors of a particular form are relevant to 
obtain the optimal asymptotics. 

Proposition 11 Let W : R + — > R + be a radially symmetric continuous and 
a.e. differentiable function. Assume that a.e. r G [0, R), 

-r < rW{r) < 3W + 2r . (20) 

Then, for all ip G H%(B R , C 2 ), 

f T f T 
/ 7 7777 W ' ^V| 2 dx > I —- \d r ip\ 2 dx , 

and the optimizers are radially symmetric and of the form ip — ( v ff ^ . In 
particular, ifW(0) > and W is nondecreasing near 0, (20) holds true for R 
sufficiently small. 

Proof. Let r = \x\. By d r , we mean V • -. For all 2-spinor ip with compact 
support in the ball B R , using (cr • ^) 2 = 1, we have 



/ 

J! 



R3 ( r + W) 



|cr • Vy?| 2 dx 



= [ 

Jm 



R3 ( r + W) 
r 

3 (r + W) 
r 



x 



{cr--){(T-V)ip 
d r ip cr ■ Lip 



dx 



dx 



I ( (\dM 2 + \\(T ■ L(p\ 2 ) dx 



~ / 7 — -T77T dr([ < ip, (T ■ Lip >] dr 

Jo (r + W) \Js* J 



= r 

Jm 



dx 



r + W) 

«3 {r + W) (lW 2 + ^k-^l 



Now, if we choose ip belonging to the class of spinors generated by the eigen- 
f unctions of cr • L with eigenvalue n, we notice that 
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,2 „ / „2 



|^ T .V?| 2 (ix + I n + ~|( — -777) ) M 2 <ir 



" L (r + W) Wl + { " r (r + " " J ^ ' 

which implies that for all ip supported in B R , 

f , ~~ . \a- Vip\ 2 dx > f , r Tjr . \d r <p\ 2 dx , (21) 

and the optimizers for this inequality are radially symmetric. 

Indeed, remember that the spectrum of (1 + a • L), is the set {±1, ±2, . . . }. 
Hence, n G {. . . , —3, —2, 0, 1, 2, ... } But our assumptions imply that the min- 
imum of (n 2 + 2n)(r + W) — (1 + W')nr on B R is nonnegative for n 7^ 
and for n = 0. Hence, the optimizers for (21) correspond to spinors which 
are eigenfunctions of 1 + a ■ L with eigenvalue 1 (n = 0). These spinors are 
radially symmetric and their second component is equal to (see [19]). 

The last assertion of the proposition trivially follows from the fact that W 
having a finite limit at 0, lim r _ >0 + rW'(r) must be equal to 0. □ 



Step 2 : We prove a relation between Hardy-like inequalities for the Laplacian 
and for the Dirac operator in the radially symmetric case. 

Consider a function W : K + — > 1R + such that W/r 3 is integrable at infinity 
and define a new variable 



1 r 2 

J^°°(s + W) s- 3 ds ~ ff°°(tr + W(tr))t- 3 dt 



Now, for any u G C^°(1R + ,M), we define q(y) := w(r), where y and r are 
related by the above change of variables. Then straightforward computations 
show that the inequalities 

r+00 ^3 r+00 

/ — \u\ 2 dr> / Wr\u\ 2 dr (23) 

Jo r + W ~ Jo 



and 



/ y 2 W\ 2 dy> V\q\ 2 dy (24) 

JO JO 

are equivalent, with V given in terms of r = r(y) by 

V(y)= 2 Z^Y n = 7 ^ f / +0 °( tr + ^( tr )) r3 ^ 2 • ( 25 ) 

v / y 2 (H^(r)+r) (r + W(r)) Ui J 
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Proposition 12 Let W : M + — > R + &e snc/i that W/r 3 is integrable at infin- 
ity. Then Hardy-like inequalities (23) and (24) are equivalent, with W and V 
related by (22) and (25). 

Remark Note that when dealing with functions which are compactly supported 
in a fixed ball, the behavior of W near infinity is irrelevant, since W can be 
modified outside the ball, without changing the integrals in the above inequali- 
ties. In particular, this is the case when searching for the optimal asymptotics 
near the origin of the functions W for which (9) holds. 

Step 3 : Let us focus now on improved Hardy inequalities for the Laplacian. 
Compared with (23), Inequality (24) is easier to deal with, because the po- 
tential appears only in the r.h.s. In [5,6,11] (see also [1,2,3,4,18] for related 
results) we find the following optimality result : 

Theorem 13 [11] The optimal asymptotical behavior near the origin for po- 
tentials V for which the Hardy-like inequality (24) holds for all q G C^°(1R 3 , R) 
is given at each order by 




An elementary proof for Theorem 13 in the radially symmetric case. For com- 
pletion, let us give a simple proof of this result. This can be done by using the 
same kind of changes of variables as those used in the proof of Lemma 9. 

Let a > 1 be the constant which appears in the definition of Xi and take 
R < e a . For all u e Hq(B r ), for every k > 1, define the functions g k by 

9k(s) := \fsg k+1 {t) , t = -^X/J) 
A simple computation shows that 

pR | rR r+oo 

/ r 2 \u'\ 2 dr = - \u\ 2 dr + \g\ \ 2 dy . 

Jo 4 Jo Jx^iR) 

With the notation t = t(s) = l/X^l/s) = a + logs, it is clear that s|| = 1. 
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From the definition of g k +i, we get, for any k > 1, 
Moreover, for any A > 0, 



2^ »»+■(*) + v^|^ + iW- 



#fc+i(*(s)) 



v^^+i 



ds 

2 



r+oo 2 



r+oo 

r+oo 
IX- 1 (A) 



g' k+ i(t) 



dt 



+ O0 



g'k+iit) 



Taking A > small enough, this means 

/; M i*wi 2 *=i/;; A) k +1 wi 2 *, 7xr , (/4) 

since gk+i has a compact support in (0, +oo). Thus 
where by we denote the inverse function of 



it 



-1/2 

rXi(r) • • -X k (r) ) g k +i 



and 



(27) 



/\ 2 |«f dr = \ l R (l+X* + ~ - + XlXl ■ ■ -Xl) \u\ 2 dr+ f + °° \g' k+l \ 2 dt . 



The asymptotical optimality shared at every order by the functions defined 
in (26) follows from the fact that for every A > 0, 



inf 



r+oo 

Ja 


g' 


2 dt 


r+oo 

Ja 


g 


2 dt 



= . 



Hence, there exists functions u such that the first term in the r.h.s. of (27) is 
negligible w.r.t. the second one. □ 

Corollary 14 Let W : M + — > R + be a radially symmetric continuous and a.e. 
differentiable function satisfying (20). Then, the optimal asymptotic growth at 
the origin for all functions W for which (9) holds in Hq(Br, C 2 ) is that of the 
function 

1 /+oo \ 

Woo(s) = l+g(E^i 2 (*)---^(*))- (28) 
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Proof. If W violates the asymptotics given by (28), a tedious calculation us- 
ing (25) shows that the corresponding potential V violates the optimal asymp- 
totics given by (26). □ 



3. 3 Optimal functions 



The first part of Theorem 4 is proved by Lemma 7, Proposition 11 and Corol- 
lary 14. For the second part, we have to match optimal functions near the 
origin and near infinity. 



+oo 



1) According to Corollaries 10 and 14, for n:— \ ^Xl(s) • • -Xj(s), 

j'=i 

1 +oo 

W := mm{W + ' 1+n , W^' 1 "™) = 1 + - ]T Xf(s) ■ ■ ■ X](s) + o(s) as s -> + 

8 3=1 

is optimal near the origin (see Appendix A for more details). A simple com- 
putation shows that W becomes smaller than 1 for any s > R, for some 
RE (0, 1). A first example of a function W > 1 which has optimal behavior 
near the origin is therefore given by W\ := max(W, 1) = min(W^ + ' m+ , W~ ,m ~), 
with 

1 ± n if s < R± , 

m± = < 

1 if s > R± , 

where [0, R±] is the support of W ±,1±n — 1. 

2) On the other hand, if we compute 

W := min(W + ' 1+fl , W^-' 1 "™) with n := ^- - 1 , 

4s 

we notice that W > 1 for all s > T for T = ^[(4096-i92\^i7) 1/a +4(4+(64+3%/4i7) 1/3 )]. 
Numerically, one finds T ps 0.866876... Hence, W2 := max(l,VF) > 1 is an 
example of a function > 1 which has an optimal behavior at infinity : 
W(s) ~ s as s -> +00. Note that = min(H/ + '™+, *-), with 



1 if s < T± , 

1 ± n if s > T± , 



where [T ± , +00) is the support of W ± ' 1±n — 1. 

The function W2 has an additional nice property : since for s large, W2 ~ s +g~) 
if we scale Inequality (9) keeping the L 2 -norm constant, on one end of the 
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scale we obtain Inequality (4), while on the other end we find the uncertainty 
principle / classical Hardy inequality (1). 



3) Now we prove that one can optimize the behavior of W near and near 
infinity simultaneously, with W > 1 on (0, +oo). 

Case + : We take a large enough so that the function W + := max(W / 00 , W +,1+n ) 
is well defined, continuous in (0, +oo) and satisfies : 

W + = in [0, R] , 

W + = w + ' 1+h in [R, +oo) , 
W + (R) > 1 , 

for some R > (numerically, a > 5 is enough). This amounts to define W + 
as W + ' m +, with 



ra+(sj 



m if s < R , 

W +,l+n j f g > R ? 

where m is the solution of the O.D.E. problem 

sm' = -2m + s (m 2 - 1) + W<x>(l + , m (0) = 1 • 
The existence of m is proved in Appendix B. 

Case — : This is dealt with in the same manner, by patching this time 
and W~' x ~ n in an appropriate way. 

The function W3 := min(W^~, W + ) satisfies all the properties stated in Theo- 
rem 4. 

In all the above examples (where W e£ 1) the functions m± have discontinuities 
and C(R) < 0. Indeed, this has to be the case whenever W 1. Let m± be 
defined by (15). According to Proposition 8, m± > 1 in a neighbourhood of 
s = + . Using W > 1, we get 

sm' ± > (m — l) 2 , 

and an easy O.D.E. argument shows that m± cannot be globally defined, so it 
must have a discontinuity. The arguments used in the proof of Proposition 6 
allow us to conclude. □ 
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Appendix A : Properties of the functions X k 



Let a > 1. Define Xi(s) := a _l for any s G (0, e a 1 ), and, by induction for 
any k > 1, X fe +i(s) := X^X^s)). Note that 

< s < e a - x =^ < Xi(s) < 1 < e 0-1 , 

which implies that s*(a) = liim _ >+oo Xfc(s) G (0,1) is independent of s (the 
limit is unique since d 2 X 1 /ds 2 changes sign only once on (0,e a )). Then 

V 2/ \ i v -i dX k+1 1 dX k 

s " Xl(s) and sX ^— = Xk+l{s) - sXk l^- 



Let vr fc (s) := n X i( s ) and a k(s) := J2 n i( s )- Since sX h+i^%T = n k+i, it 
follows that s ^7 = 7T fc cr fc . By definition of X fc , X fc+ i(s) = -^(^^x^s) and 



<J k+1 (s)=t(l+a k (t)) a nd S ^±i= ft ^(t) + a k (t) + lV 

/|t=Xi(s) as \ at J \t=xi(s) 



Using the two above identities, we can prove by induction the following for- 
mula : 

Lemma 15 For any k > 1, for any s G (0, e a ~ l ), 

2s^(s)-al(s) = ^(s). 



We may now pass to the limit k — > +oo. Let a(s) : = Dj=i n j( s ) '■ 

as j=1 

With the notations := 1 + § ttJ(s) and n(s) := \ a = \ ^(s), 

this means 



Corollary 16 For any s G (0,e a 1 ) 7 
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Appendix B : Solving a singular O.D.E. 



Here we solve the differential equation 

sn' — n 2 — 2 sn — sn 2 . . _ . . 

W OQ -l = , n = 1 , 29 

oo 

in an interval [0, 5], 5 > 0, small, with W^s) = 1 + - ^ x i( s ) 2 • ' -X k {s) 2 . 

8 k=i 



Proposition 17 There exists 5 > such that (29) has a continuous solution 
in [0,5). 

Note that this problem is a limiting one in the sense that there is no function W 
more singular than at the origin, for which the above problem can be 
solved with continuity at the origin. 

Proof. Let C, 5 be two positive constants and define the set 

X c>5 := {u G C([0,5]) : limsup| M (s) log(s)| < C} . 

Let us write n :— n (1 + w) — \ Y^=\ -^i( s ) • • • Xj(s). Then n is a solution 
to (29) if and only if w is a solution to 

w = fo + fiw + f 2 w 2 , (30) 

where fo, fi, ji have the following behavior near : 



1 

4 s | logs| 2 
1 

s I logs| log(| logsj) 

1 

1 s I logs] 

In order to solve equation (30) together with the initial value w(0) = 0, we 
introduce the map T : Xc,s — > Xc,s defined by 

Tw(s) := f\f + fiw + f 2 w 2 )dy, 
Jo 

and look for a fiwed point. By choosing C > 1/4 and 5 < 1 small enough, 
T maps X(7 j( 5 into itself and it is a contraction. So, there is a unique solution 
of (30) in Xc,s which means that (29) has a unique continuous solution n in 
the interval [0, 5]), with n/n — 1 in X c> 8, such that n(0) = 0. 



fo = 
fi = 
/ 2 
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Appendix C : Why do we have a liminf in Theorem 4 ? 



We are going to give a qualitative example showing that only a liminf can be 
achieved. Let W := + ^ w(- — — ) where e n and s n are such that 



s n > , lim s n = , 

£ n > , £ n < +°° > 

n>0 

£ n+l *^ S n S n _)_i , 



and assume that FT is a bounded function with compact support in (0, 1). 
Then 

lim sup W(s) > 1 

s->0+ 

and the equations 

2m± ± sm' ± — sm^_ + s = W(l + m^_) 
have no solution continuous up to s = 0. 
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